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In thii3 paper the authors argue that the decision 
structures of educational institutions are mult i-- level, 
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in decision theoretic terms and that these multi-levels, 
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computer program used for the problem and the input data 
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included. (Author) ^ 



FILMED FROM BEST AVAILABLE COPY 



OO 

^.-.^ 
CO 

o 
a 




ERIC 



FORD FOUNDATION PROGRAM FOR 
RESEARCH IN UNIVERSITY ADMINISTRATION 

Office of the Vice President — Planning 
University of California 



us DEPARTMENT OF HEALTH 
EOUCATIONA WELFARE 
NATIONAL INSTITUTE OP 
EOUCATION 

Tm.s DOCUMfNT HAS BEF. N WFPWO 
Que ED EXACTLY AS 5?ECrivtD f r^CV 
THE PERSON 0» ORGANI?ATiON OCiO'N 
ATJN&iT POINTS or VI C/; OP OPi NIONS 
MATED 00 NOT NECeSSAWiLY REP«E 
SENTO*f"IClALNAT.ONALiN«=.TITuTEOF 
EDUC:*TiON POSITION OP PQUCY 



FORD GRANT NO. 680-0267A 



RESEARCH DIRECTORATE 



Charles J. Hitch 
Frederick E. Balderston 



George B. Weather sby 



President, University of California 

Professor of Business Administration 

Chairman, Center for Research in Management Science 

University of Cahfornia, Berkeley 

Academic Assistant to the President 

Associate Director, OfiBce of Analytical Studies 
University of California 



OFFICE ADDRESS 

2288 Fulton Street 
Berkeley, California 94720 
(415) 642-5490 

(List of Available Publications on Inside Back Cover) 



OPTIMALITY IN COLLEGE PLANNING 
A CONTROL THEORETIC APPROACH 

W. Gary Wagner 
George B. Weathersby 



Paper P-22 
December, 1971 



TABLE OF CONTENTS 



Page 

PREFACE ii 

LIST OF TABLES • • • • iii 

LIST OF FIGURES iv 

INTRODUCTION 1 

INSTITUTIONAL DECISION MAKING . . 4 

Analytical Description of Institutional Decision Making 11 

The President's Model: Single Decision Maker Paradigm 15 

SMALL CAMPUS PLANNING MODEL ..... 19 

SYSTEM DYNAMICS OF SCPM . 20 

Faculty . 20 

Students 21 

Space 22 

Money . , 22 

CRITERION FUNCTION . . . 25 

Student/ Faculty Ratio 25 

Faculty Mix 26 

Space 27 

Student Mix 27 

Monetary Balance 28 

Terminal Conditions 29 

SUMMARY OF RESULTS ...... 30 

Comparison of Optimal versus Sub-Optimal Results 34 

Comparison of Results from Various Policy Alternatives 40 

Conclusion 46 

APPENDIX .47 

Solution Procedure ... 47 

Computer Program . 48 

Data Inputs for Multiple Runs 60 

Form of Input Data 61 

Data Structure for Multiple Runs 62 

BIBLIOGRAPHY 63 



ERIC 



1 



PREFACE 



This is one of a continuing series of reports of the Ford Foundation 
sponsored Research Program in University Administration at the University 
of California, Berkeley. The guiding purpose of this Program is to under- 
take quantitative research which will assist university administrators 
and other individuals seriously concerned with the management of univer- 
sity systems both to understand the basic functions of their complex 
systems and to utilize effectively the tools of modern management in the 
allocation of educational resources. 

In this paper we argue that the decision structures of educat5onal 
institutions are multi-level, multi-decision-maker hierarchies which can 
be described and analyzed in decision theoretic terms. Furthermore, we 
show that these multi-level, multi-decision-maker hierarchies can be re- 

t duced to equivalent one-level, one-decision-maker formulations, which 

\ 

V^an be solved either analytically or numerically by the techniques pre- . 
sdnted in this paper. An illustrative example is given which first 
identifies and then solves for a set of optimal resource allocation and 
policy .^de visions. A listing of the computer program used in this problem 
and the input data specifications are included in the Appendix. 
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INTRODUCTION 

Institutions of higher education currently face a number of major 
policy choices which will largely determine their character for the next 
twenty-five years. The tremendous expansion of American higher education 
in the last twenty-five years was driven by burgeoning enrollment growth 
and by massive federal commitment to doctoral production in the sciences 
and technologies. Both of these forces are abating rapidly. Nationally, 
enrollments in higher education are forecasted to peak in 1980, then de- 
cline until the late 1980' s and not approach the 1980 level until after 
1995. Many schools are now experiencing enrollment levels below their 
previous expectations. This is not a short-run phenomenon; rather, cur- 
rent enrollment shortfalls are harbingers of the next twenty-five years. 
Colleges and universities must learn to survive and to prosper with a 
decreasing demand for their services. 

It is far less likely that in the future the federal government will 
rescue the expectations of higher education as they have done in the post- 
Sputnik era. The United States will probably have a surplus of highly 
trained scientific and technical manpower for at least the next decade 
without major additional federal expenditures [Erode (1971)]. The re- 
duced rate of undergraduate enrollment expansion will drastically reduce 
the number of new teachers needed in colleges and universities, thereby 
reducing the future demand for additional Ph.D.'s and, therefore, the 
need for large doctoral programs. Furthermore, the federal priorities 
have shifted from scientific manpower to equality of student access and 
the quality of educational experience. Both of these major federal 
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objectives impact undergraduate education far more than graduate programs 
and they move counter to many institutions' prestige and elitist orienta- 
tions . 

' Ivhile the demands for educational services by students and govern- 
ments will probably be decreasing in the coming decades (first in rate of 
growth and then in absolute number), the costs of educational institutions 
continue to rise. As Cheit [1970] has pointed out, a significant number 
of America's colleges and universities are headed for financial difficul- 
ties and current institutional rigidities preclude those cost adjust- 
ments necessary to maintain fiscal viability-"^ Furthermore, the techno- 
logy of education has changed very little in the last three or four 
decades; indeed, some would argue that educational technology has changed 
very little since Socrates. In essence, there has been no observable 
productivity increase in American post-secondary education in the last 
four decades [O'Neill (1971)]. 

If it were not so painful, we might examine with considerable in- 
tellectual interest the experience of public and private- eleemosynary 

institutions beset by diminishing demands for services and risiag costs, 

2 

increasing institutional rigidities and no productivity increaso.s. 
Unpleasant as it may be, educational administrators are having to ask the 



"^The "varieties of the financial crisis" are explored by Balderston 
[1971] in a recent paper prepared for the American Council on Education. 
2 

This description closely resembles the experience of the American 
railroads. One inciteful observation on the decline of the railroad 
companies was that unfortunately railroad managers viewed their industry 
as "railroads'* versus "transportation." At the time of burgeoning new 
modes of transportation, the railroad companies were in an excellent 
position to diversify and expand—but that was neither their tradition nor 
their self-concept. Are our schools in the "formal instruction" industry 
or the "education" industry? 



tough questions: What are our objectives? How would we know if we 
achieved them? How can we reallocate resources to be more productive? 
What activities are really essential for an educational institution? 
Who should make these decisions; and many others? 

The purpose of this paper is to look at institutional resource 
management decisions in the context of institutional goals and objec- 
tives. After describing one view of institutional decision making, we 
present a simple yet comprehensive mathematical model which explicates " 
the interrelationships of major institutional variables. Sample data 
are then used to derive resource allocations which would be optimal for 
the institution. The use of this model in educational policy analysis 
is then discussed before presenting our conclusions. 



INSTITUTIONAL DECISION MAKING 

The decision making structures of educational institutions are as 
diverse as the institutions themselves. Some colleges and universities 
are highly authoritarian while others are highly egalitarian; some insti- 
tutions are ruled by presidents and others by committees • Some educa- 
tional systems have many layers of administrative superstructure while 
others do not. There is often little resemblance between the organiza- 
tional structure and the decision or power structure of an institution* 
Often individuals vith no delegated authority have great influence on 
decisions. 

While these complex interactions have been analyzed from many per- 
3 

spectives, we have chosen to analyze institutional decision making from 
a decision theoretic basis. Initially we distinguish between the values 
used in arriving at policy decisions and the authority structure in 
which the decisions are made. Focusing first upon the structure of de- 
cisions, we observe that the decision structures of most educational 
institutions are hierarchical, with students, individual faculty members, 
department chairmen, deans, provosts, and presidents playing different, 
but important, decision making roles* 

These roles are distinguished primarily by the variables each level 
can control. For example, students decide which of the available courses 
they will take; faculty decide how to allocate their time between formal 

3 

Wildavsky [1964] looks at the resource allocation process in govern- 
ment from a political theory perspective; Downs [1967] and Braybrooke and 
Lindblom [1963] view bureaucratic decision making as a behavioral and organi 
zational process; Glenny [1969] and Palola [1970] approach educational 
decision making from the perspective of governance. 
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instruction, preparation, informal meetings with students, research, 
committees, community service, professional advancement, and leisure; 
department chairmen decide, with consultation to be sure, the course and 
committee assignments of faculty, the allocation of support services, 
recommendations on salaries and promotions, and curriculum proposals; 
deans allocate new faculty positions to departments, increasingly will 
reallocate faculty positions between departments, determine salaries and 
promotions^ establish departmental budgets, endorse curriculum changes, 
and approve research programs; provosts or presidents in turn allocate 
faculty positions and budgets between schools, review or approve salaries 
and promotions, recruit deans, approve curricula and academic programs. 
Table 1 summarizes some of these distinctions. 

Another characteristic of the hierarchical structure of educational 
decision making is the direct interrelationships of the various decision 
making levels. As illustrated in Table 1, the control variables at one 
level often become constraints at the next lower level. For example, 
the president can allocate faculty positions to the various schools in 
his institution to the limit of his budget. In turn, deans can allocate 
faculty to departments up to the limit permitted by the president's bud- 
get. What was initially a decision to the president later becomes a con- 
straint to the dean. 

Another component of our analysis of institutional decision making 
is the distinction between the implementation structure and the decision 
structure. The implementation structure is usually reflected in the 
institution's organization chart; it is the array of deans, department 
chairmen, accounting officers, purchasing agents, budget officers, ad- 
missions officers, registrars, librarians, and all the other functional 
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TABLE 1 

Examples of Interrelationships of Institutional Decision Makers 



Decision Makers 


Control Variables 


Constraints 


President/Provost 


Budgets 

Faculty Positions 
Program Approval 


Income or Appropriations 


Dean 


Departmental Budgets 
Faculty Positions 
Program Approval 


President's Budget 
President's Budget 
President's Approval 


Department Chairmen 


Faculty Assignments 
Support Services 
Salaries and Promotions 
Curriculum 


Dean*s Budget 
Dean's Budget 
Dean's Budget 
Dean's Approval 
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specialists who keep an institution running effectively. On the other 
hand, the decision structure is rarely reflected in a school's organiza- 
tion chart. At issue here is who is responsible for what decisions and 
how are the recommendations for these decisions made. 

The operations cycle of an institution is ilJustrated in Figure 1. 
Once made, a decision is communicated to the implementation structure 
where functional specialists establish the operating policies and pro- 
cedures which actually move the organization in the desired direction. 
These implementation managers need operating data for their effective 
functioning. For example, the accounting officer needs payroll infor- 
mation to process checks and charge the appropriate accounts. On the 
other side of this circle, those who are charged with recommending de- 
cisions need institutional data to evaluate past decisions and as input 
to future decisions. In addition, external data on student demand, man- 
power supply and demand, com.munity needs, attractive research areas and 
a variety of other issues are needed for decision recommendations- The 
process of synthesizing these data and institutional objectives into a 
coherent, consistent strategy for future action is nebulous if not non- 
existent at most institutions. Yet, if educational leaders are going to 
be able to deal effectively with the serious social and economic challenges 
confronting their institutions, much more attention will have to be de- 
voted to their decision structures. The approach and the mathematical 
model presented in this paper is one small step in this direction. 

In addition to the process of decision are the values upon which 
the decision is based. One of the functions of each decision maker 
is to choose the values appropriate for the decisions at hand. 
This is another way of raising the question of governance: Who will 
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FIGURE 1 



Operations Cycle of an Institution: 
Decision - Execution - Evaluation - Decision 
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decide and whose values will he use when he decides? Furthermore, what 

attributes of the educational system are important to the decision maker, 

4 

what does he consider to be the outputs of his system? Hov important 
are more undergraduates versus more graduate students, more researchers 
versus more instructors, more computing power versus more library services, 
more faculty versus more facilities, and a thousand^and-one other possible 
tradeoffs? Taken together, all of these choices and tradeoffs comprise a 
decision maker^s value system. 

These value systems also serve to connect the hierarchical decision 
systems which were discussed earlier. In many cases, the president of 
an institution is deeply concerned about the classroom environment and 
the interaction of students and faculty, even though he cannot directly 
control any of the operative variables. However, the president often 
makes his budgetary and faculty allocations with their educational con- 
sequences in mind and adjusts his allocations to correspond to his assess- 
ment of the educational use to which these resources are put. In other 
words, the consequences of decisions at a lower level are important to 
decision makers at higher levels. 

There is a circular flow of information in a hierarchical decision 
system: decisions are passed downward and value signals are passed upv;ard . 
These in turn affect the decisions which are passed downward in a subse- 
quent cycle, as shown in Figure 2. It is this two-way flow of information 
that makes delegated authority operational and renders a decentralized or 
hierarchical system controllable; this notion of decentralized control will 
be explored in more detail shortly. 

4 

The outputs of higher education have received increasing attention in 
the last few years; see Lawrence, Weathersby and Patterson [1970], 
Breneman and Weathersby [1970] and Huff [1971]. 
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Constraints 
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Values 


Consequences 


Decisions 


Environment 


Constraints 









Department 
Chairmen 




FIGURE 2 

Interrelationship of Decisions and Values 
in a Hierarchical Decision Structure 
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Analytical Description of Institutional Decision Making 

This conceptual analysis of educational decision making can be made more 

precise by describing the decision interrelationships in mathematical terms. 

To begin with, we need some definitions:^ 

u.(t) == the v6iCtor of decision variables available to decision 
1 

maker (DM) i in period t ; 

z^(t) = the vector of predetermined variables impinging upon the 
system relevant to DM i in period t ; 

x^(t) = the vector of consequences (or state variables) in 

period t resulting from the decisions of DM i and rele- 
vant exogenous influences. 

The relationship of consequences to decisions (or output to input) is often 

called the production function: 

x^(t+l) = f^(x^(t), u^(t), 2^(t), t). (1) 

Finally, the value to DM^ of making a decision ^j^(^) when confronted 
with the predetermined variables ^j^^^) written 

V.(x.(t), u.(t), 2.(t), t). (2) 

Expressions (1) and (2) describe the horizontal flows £;hown in Figure 2 
at each decision making level. 

The decision problem faced by each administrator is to maximize his 
own values subject to his constraints of authority and resources and sub-^ 
ject to the responsiveness of his system to the application of policy or 
resource decisions. Furthermore, a decision maker often looks several 

"^For a complete exposition of this approach and the motivation for 
these definitions,, see Weathersby [1969a, 1969b]. 
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years in advance and wants to Tnaxiinize his values ovei a planning horizon 
of N periods. We may write this dec:*=ion problem as:^ 



N 

max {j - I V. 

u^(o), . . • ,u^(N-l) 1 t=o 



x^(t), u^(t), z^(t), t 



subject to: 



(3) 



x.(t+l) = f. 
1 1 



X^(t), U^(t), 2^(t), t 



t=o ,1. . . ,n-l 



(1) 



X . (o ) known 
1 



and 



x.(t) 
u.(t) 



< b.(t) 



t-o ,1. . . ,n-l 



(4) 



In' this formulation, C is the constraint function and ^-j^(t) 
resource and other constraints relevant to DM. in period t . 

In general, there is a solution u^(t), t=0,l, . . . ,n-l which maxi- 
mizes the overall value function, J , provided the necessary and suffi- 
cient conditions are satisfied.^ Furthermore, an optimal solution is in 

g 

general a function of all preceding variables. 

u^(t) = g^ u^(o),... ,u^(t-l), x^(o) , . . . ,x^(t) , z^(o),...x^(t), 
b^(o),*..b^(t), tj . (5) 



can 



The general form of the N period value function is V^^x^(o), 
x.(l),...,x.(N), u (o), u.(l),...,u.(N-l), z.(o),...,z.rN-l)l which 

be separated into V^^x^(t), u^(t), z^(t), tj + V^^x^(N), Nj by the 

assumption of weak separability (Weathersby [1969a]). N->tice that present 
value discounting is a special case of the time functiou of V(x,u,z,t). 

'^The solution procedure will be described in more detail later and 

the algorithm used in this study is described in the Appendix, 
g 

See Aoki [1969] for a discussion of general recursive solution. 
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The consequences which result from an optimal decision sequence can be 
Calculated from equation (1); 



x*(t+l) = f.lx.(t), u*(t), z.(t), tj 



= h* 

1 



X.. (o) , . . . ,x^(t) , z^(o) . . ,z^(t) , b^(o) . . ,b^(t) ,t 



(6) 



after equation (5) is substituted for u^t) . In other words, when the 
parameters of the decision problem are known, i.e., equations (l) , (3), 
(4), and (5), one can replace the decision problem by equation (6) which 
describes the consequer -^es cf an optimally controlled system. 

In a strict hie chicai decision structure, there is one such deci- 
sion problem for each dicisi.-u tiaker. In the three level hierarchy shown 
in Figure 3, n deans /eport to the president and m^ department chair- 
men report to the i^^ dean. The interrelationships are: 

1. At levels 1 and 2, the predetermined variables 2(t) and 
resourcf constraints b (t) can be controlled or influenced 
at the lext higher level. For example, 

^,1^'-'^^ = ^l,lK,l^^)' ^1.1^'^' "2.1^'^' ^2,1^')' 

2. The values associated with the consequences of decisions by 
decision-makers at levels 2 and 3 can include the decisions 
and consequences of lower level decision-makers. For example: 

^,1= ^3,1^')' ^3,1^'^' ^2,1^^)' x,^,(t),....x^^^,(t). 

S^l(t),...,U2^^(t). t 

The first interrelationship describes the downward flow of decisions either 
directly or indirectly while the second form of interrelationship describes 
the upward flow of accountability or value. Both of these interrelationships 
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level 3 
(President) 



level 2 
(Deans) 



level 1 

(Depart. 

Chairmen) 




FIGURE 3 

Array of Decision Makers in Three Level Hierarchy 



are necessary for a controllable hierarchical system; however, they are not 

sufficient conditions for total systems controllability. 

The basic strategy of solution of a strict hierarchical decision 

structure (i.e., interrelationships only between adjoining levels) is to 

reduce the structure down to a one decision maker problem by folding up 
9 

from the bottom. This approach would replace each decision making node 
at level 1 in Figure 3 by his corresponding optimal decision function 
g . , and the corresponding optimal production function h. . . . For 
each level 2 decision maker, the m^ departmental chairmen's decision 



For a complete presentation of this approach, see Weathersby [1969b] 
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problems are replaced by 2Tn^ vector equations for the g and, h func- 
tions. These 2in^ vector equations then are effectively production functions 
' t h 

to the i dean which augment his own production function. Now wo can fold 
the dean's decision problems up to the presidents' level by the same tech- 
nique. Thus, we can collapse a multi-layer, multi-decision maker hierarchy 
to an equivalent one-decision-maker problem. Correspondingly, if we can 
solve the one-decision-maker decision problem, we can conceptually solve 
the multi-layer, multi-decision maker problem. Therefore, the remainder 
of this study will focus on the single decision maker problem. 

The President's Model: 

Single Decision Maker Paradigm 

The basic decision problem of the president is to maximize the achieve- 
ment of his own values, or the values he chooses to operate with as presi- 
dent, subject to resource limitations and the responsiveness of his insti- 
tution. The formal statement of this problem was given previously in 
equations (3), (1), and (4). The three major components of the problem 
are: (1) the president's value function; (2) the institutional response 
or production function; and (3) the resource, legal and other constraints. 

One obvious difficulty with the decision theoretic formulation is 
that generally presidents, and other administrators, cannot articulate 
their value function* Mofit of us are not trained to think in terms of 
multi-attributed utility functions and, therefore, any approach which re- 
quires a mathematical description of a decision maker's value function is 
destined to grave difficulty if not failure. 

There have been two major techniques for circumventing the assess- 
ment difficulties associated with a full description of the value function. 
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Geoffrion and Dyer [1970, 1971] have shown that one need only assess the 
local gradient of the value function at the current operating point. In 
their work, they ask a dean or department chairman to select one of his 
variables as a numeraire and then assess the pairwise tradeoffs of all the 
other policy relevant variables with respect to the chosen numeraire. Tlais 
is the local gradient which shows an improving direction along which the de- 
cision maker selects a new and improved operating point. At this new point, 
however, the local gradient must be assessed again because it is generally 
different at eyery point on the utility surface. In other words, the 
Geoffrion and Dyer approach replaced a global assessment of the multi- 
dimensional value surface with a series of local assessments of the tangent 
plane, which is a much easier task. This requires interaction between the 
decision maker and the mathematical programming algorithm because the path 
along the value surface is unpredictable a^ priori . 

A second approach to the reduction of the dimensionality of the value 
assessments is to express the decision-maker's objectives in terms of 
targets* This is the approach used in the study reported in this paper. 
For simplicity of exposition, consider a president's value function that 
is defined over only the consequences of state vector x , i.e., V(x) , and 

is 

that the president wants to achieve a most desirable level of x , say x 
In other words, the president believes that the optimal state of his institu- 
tion would be a student enrollment of 10,000 with 1,000 faculty members, 
600 of whom would be tenured, and so forth. 

The key to the target approach is that if the institution is initially 
reasonably near the desired targets, then the general utility maximization 
problem can be expressed by an approximately equivalent loss minimization 
^^jblem where the loss function is quadratic, independent of the form of 
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the utility function as long as the utility function is twice differen- 
tiable, i.e., smooth and continuous. Although we may know nothing more 
about the president's utility function than that it ought to be concave and 
smooth, we do know that in the neighborhood of his targets his loss 
function is quadratic to second order* 

By choosing targets x , a decision maker indicates that 

V(x*) > V(x) for all x . (7) 

For X near x > we can expand V(x) about x by Taylor series 

V(x) = V(x*) + wl (x-x*) + l/2(x-x*)^ V^V j (x-::'^) 

X X 

Higher Order Terms. (8) 
If V(x ) is a maximum, as indicated by expression (7), then the local 

•k 

gradient must be zero at x and the second right hand side term in equa- 
tion (8) must be zero. Furthermore, the second derivative of V must be 
negative definitive for x* to be a strict maximal point. Therefore, the 
third right hand side term in equation (8) must be negative for all x . 
This argument proves that, to second order, 

max V(x) = min 1/2 Ux-x*) ^V^v| (x-x*)} . (9) 
X x ^ x* ' 

One point of indeterminancy remains in equation (9) ; in general, the 
matrix of second partial derivatives of V is not known. Two approaches 
may be used here. One can ask the decision maker to choose a numeraire 
and assess the relative pairwise comparison losses that he would exper- 



^ we can assess the relative relationship between the first and 

pn¥/^"ond derivatives of a decision^maker ' s utility function by a discussion 
^jj^^his risk aversion (Pratt [1964]). 
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lence at x and use this one set of assessments in place of V V . Alter- 

2 

natively, one could recognize that the magnitudes of V V change the rela- 
tive shape of the quadratic loss structure but not its minimuTn, which is 

X . Near the minimum, the solution to (9) is often insenritive to the 

2 

global shape of V V and a much simpler procedure is possible: namely 

2 

choose an arbitrary weighting matrix K in place of V V such that the 
magnitudes of loss of one unit variation in every dimension are identical. 
Both of these approaches require minimal assessment » 

In summary, we have argued that the decision structures of educational 
institutions are multi-level, mult i-decis ion-maker hierarchies which can 
be described and analyzed in decision theoretic terms. Furthermore, we 
have argued that these multi-level, multi-decision-maker hierarchies can 
be reduced to equivalent one level, one decision-maker formulations. In 
turn 5 these single decision-maker problems can be solved either analyti- 
cally or numerically by the techniques discussed in this and the following 
section. We now proceed with the formulation and solution of a specific 
decision model and discuss its implications* 
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SMALL CAMPUS PLANNING MODEL 

The concepts of the previous section are illustrated in this section in 
a specific analytical modeling context. For the purposes of exposition, we 
have focused on the instructional program of an institution partly because 
this seems to be an area of great interest to most colleges and partly be- 
cause instructional activities have far more in common among institutions 
than the various research and public service programs. The paradigm of this 
model is the liberal arts undergraduate institution or that component of a 
major university. 

The Small Campus Planning Model (SCPM) is designed to provide a con- 
trol theoretic solution to the problem of finding an optimal sequence of new 
student admissions, new faculty hires, and new physical construction over 
an N-year planning horizon. It assumes that the flows of students, faculty, 
construction, and money can be characterized by linear dynamic equations 
and that the campus administrator's preferences for student and faculty mix, 
for space, and for solvency are sufficiently close to the ins ^.i tut ion' s cur- 
rent experience that actual deviations from targets can be adequately ex- 
pressed in terms of quadratic penalty functions. 

The model is still in the investigatory stages and wil] undoubtedly 
undergo further revision before it is considered a finished product. Ul- 
timately it is hoped that SCPM, because of its minimal data requirements 
and low implementation and calculation costs, may serve as a useful plan- 
ning device for college administrators who have neither the funds nor the 
data base to support implementation of other, more complex, models. "^"^ 
^ "^"^For a structural comparison of other recent analytical models for 




'versity planning, see Weathersby and Weinstein [1970]. 
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SYSTEM DYNAMICS OF SCPM 



The generalized form of SCPM's dynamics may be characterized by 



x(t+l) = X (t) + u(t) + 2(t) 



(1) 



where : 



x(t) 



n-vector of state variables at time t 



u(t) 
z(t) 



m-vector of control variables at time t 



n-vector of predetermined variables at time t 




n X n matrix of transition coefficients for period t 
n X m matrix of transition coefficients for period t 
n X n matrix of transition coefficients for period t 



It is assumed in SCPM that = F , = G , 



and ^ H for all 



t , i.e*, that the transition matrices are not time dependent. This is 
not necessary for solution, but facilitates estimation and reduces sub-- 
stantially the data requirements. Furthermore, ^ hese matrices are not 
Markovian, i.e., the row sums do not total 1.0, because the absorption 
states for students, faculty, space and money are excluded. Under the 
stationarity assumption, equation (1) becomes 



The state, control, and predetermined variables as defined in SCPM are 
divided naturally into four groups: students, faculty, space, and money. 
We consider each separately for purposes of exposition. 



We define the variables ^]^(^) ^4^^^ number of full- 

time equivalent full professors, associate professors, assistant professors 



x(t+l) = F X (t) + G u(t) + H z(t) 



(2) 



Facul ty 
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and instructors who were in the institution last year (at time t-1) and 
who remained in the system at the start of this year (period t) . Similarly, 
we define the variables ^]^(^) ^4^^^ number of faculty who 

are hired at corresponding raaks at the start of period t » Then, ex- 
pressing the matrices F and G by their elements f,^^ and g^^ , we have 

x^(t+l) = + ^12^2^^^ ^ll^'l^^^ ^ ^12''2^^^ 



^4(^+1) - ^44^4^'^ ^44^^'^ ' 

Here, f_^^ is the promotion rate of faculty from level j in period t 
to level i in period t+1 , f.. is the continuation rate for faculty 
at the same rank, and g_ is the promotion and continuation rates of 
new faculty who were hired in period t at level j and who are at level 
i in period t+1 . 



Students 

We define the variables Xj-(t) to Xo(t) to be the number of continu- 

ing freshman, sophomore, junior and senior students at the start of period 

t . Similarly, we define the variables u-(t) to Uo(t) to be the number 

J o 

of students admitted to the corresponding student levels at the start of 
period t . Once again, we may write the scalar equations 

x^(t+l) = f33X^(t) +g55U3(t) 

xg(t+l) = fg^x/t) + f33X3(t ) + ggyu/t) + g88^8^'> " ^^^^ 

The coefficients f.. and c.. have the same advancement and retention 
interpretations as before. Attrition of faculty and students is accounted 
^ for by omission of state variables corresponding to the **out" state. For 
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manpower planning or other purposes one could define two additional states 
of successful degree completion and "stopping out.** This would provide 
Specific degree output information and render the student system Markovian. 

Space 

We assume that physical construction takes an average of four 
years to complete once it has begun. A conscious simplification at this 
stage is the assumption of fully interchangeable space types and uses. 
Equation (5) could be repeated for each space type if the additional de- 
tail would be worth the additional cost. SCPM also assumes a constant 
depreciation rate of a = l^~^i2 12-' * Accordingly we define 

Ug(t) = amount of new construction measured in Assignable 
Square Feet (ASF) which begins in period t ^y- 

Xg(t) = ASF begun in period t-1 

x^^Ct) = ASF begun in period t-2 

x^^(t) = ASF begun in period t-3 

x^^Ct) = ASF which is available and usable at the start of period 



Thus, 



Xg(t+1) = Ug(t) i.e., gg g= 1.0 
x^Q(t+l) = xg(t) i.e., f^Q^g = 1.0 
x^3^(t+l) = x^Q(t) i.e., f^^^^Q = 1.0 

x^2(t+l) = x^l(t) + f^2,12 ^12^^^ ^•^•'^12,11=^-^ • 



Money 

Finally, we assume that there are two kinds of funds which adequately 
describe the administrator's financial concerns: restricted funds (endow- 
ment) and unrestricted funds (operating plus capital funds). Once again. 
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these fund categories could be expanded if needed. It is further assumed 
that interest earned on endowment funds raay be allocated arbitrarily be- 
tween funds, but that the income and capital gains use policy is fixed in 

advance. If f., . is the value (cost if negative) of one unit of x. 

14, J ^ J 

and / . is the value (cost) of one unit of u, , and performing all 
1^ J J J 

calculations in constant dollars, we can define 

z^^(t) = restricted gifts in period t (estimated or assumed known) 

Zj^^(t) = unrestricted gifts in period t (estimated or assumed 
known) 

and write: 

x^3(t+l) = fi3, 13x^3(0 + h^3^^3Z^3(t) 

x^^(t+l) = I fi4^^x.(t) + I gi4^.u.(t) + h^^^^3Z^3(t) + h^4^i^z^4(t) 

While equation (6) looks quite complicated, each of its components 

is very simple and traditional. 

f^, ^ to f|, , = the average faculty salary by rank including 

direct support costs. 

■^14 5 ^14 8 ^ institutional cost per student by level 

' ' (excluding faculty salaries and direct faculty 

support costs and including student fees and 
tuition) . 

^14 9 ^14 12 ^ average cash flow cost per ASF in each 

' ' year of new construction. 

f ^ , ^ ^ = the average rate of return on endowment that 

14,13 . . 1 , 1 1- 

' IS available for operating expenses. 

f^^ = the proportion of last year's net cash balance 

' available in the current year (usually 1.0). 

The numerical values of these and the other coefficients used in the com- 
putational example are given in Table 3. Note that equation (2) is just 
the aggregate of equations (3) through (6). If the state vector x(l) is 
known and the gift funds z(t) , t = 1,...,N-1 are predicted, then the 
control sequence u(t), t = 1, . . . ,N-1 with (2) determines x(t) for all 
future periods t - i,...,N . After making various assumptions about the 
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likely levels of state aid and gifts in the future, SCPM determines optimal 
enrollment, hiring, and construction policies for a given set of institu- 
tional objectives, which are described in the next section. 

Student tuition could be included as a control variable instead of a 
predetermined factor. This would more accurately reflect the decisions of 
most private institutions and a growing number of public institutions. How- 
ever, two major problems have to be dealt with to include tuition as a uncon- 
strained Control variable. A conceptual problem is the effect of additional 

tuitr'on on student demand for attendance and on the quality , of students able 

12 

to pay the higher tuition. A minor technical problem is the non-linearity 
of the money dynamics introduced by controllable tuition. The solution al- 
gorithm given in the Appendix will accommodate both linear and nonlinear dy- 
namic systems. 

However, this does raise the issue of the validity of the linearity as- 
sumption embodied in the SCPM systems dynamics. While in any specific imple- 
mentation, the functional form of the dynamics would be an empirical ques- 
tion, there are several justifications for the use of linear dynamics in 
our example: (1) the ease of interpretation of coefficients in terms of tran- 
sition probabilities, depreciation factors, faculty salaries, etc.; (2) the 
experimental ease of formulation and modification; and (3) the lack of any 
information of a more generally useful and accurate formulation. 



See Miller [1971] for a discussion of recent attempts to estimate 
student demand functions and Jewett [1971] for a presentation of a national 
student ability - willingness to pay model and analysis. 
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CRITERION FUNCTION 



The general form of the criterion function used by SCPM is 



min 



N-1 

J = p[x(N)) + I V (x(t),u(t),z(t)) 
t=l 



(7) 



where P and are the relative quadratic loss functions darived in 

a previous section. In this study, P is the sum of four quadratic terms 
and (V^=V for all t) is a summation of nj.ne quadratic terms. For 

purposes of exposition we separate V into five sets of terms relating 
to the administrator's objectives expressed in terms of the student/ 
faculty ratio, faculty mix, space requirements, student mix, and finan- 
cial stability. F will be discussed separately. 

Student/Faculty Ratio 

One proxy measure of the amount and quality of student/faculty inter- 
action at an institution is the ratio of students to its (PTE) faculty. 
SCPM enables a campus administrator to specify a targeted ratio and then 
seeks a set of controls which minimizes the deviation of the actual student/ 
faculty ratio from the target. If we define 



r 



1 



target student/faculty ratio. 




some scalar weight. 



\ fx (t) + u (t)] , 
I fx. (O + u.(t)l , 



which is the total faculty at time t, and 



which is the total students at time t, 



i=5' 



then the first term of V may be written 



k^(r^TF^ - TS,)2 



(i) 
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If an administrator wished, this term could be expanded to reflect stu- 
dents by level and faculty by rank to describe, for example, the exposure 
of lower division students to tenured faculty. Similar expansions are 
possible in most of the terms of the criterion function but were not re- 
ported here because they do not alter the solution procedure or the basic 
utility and results of the model. 



Another proxy criterion for assessing the quality or prestige of a 
college is the mix of its faculty by rank. [See Rowe, Wagner and Weathersby 
(1970).] In the case of community colleges, or any other cases for which 
there are no ranks but rather salary schedules, we may interpret the four 
(or fewer) levels of faculty purely in terms of salary. In any case, we 
assume that the administrator has preferences over different mixes of 
faculty by level and we allow him to specify target ratios which describe 
the desired mix. 

If .we define 



Faculty Mix 




x^(t) + u^(t), i = 1,2,3,4, which is the number of FTE 
faculty at level 1 in year t , 



r 



i 



target ratios of each rank relative to the number of 
full professors, i = 2,3,4, and 



k. 



1 



scalar weights i 



= 2,3,4, 



then the next three terms of V 



are 




i = 2,3,4 



(ii) 
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Space 



Typically space needs are largely determined by either student en- 



rollment or faculty size or a combination of the two. SCPM makes the 



simplifying assumptions that space is interchangeable and available 
continuously variable amounts. For a small, homogenous college, space 
interchangeability may not be a devastating assumption because in liberal 
arts subjects rooms of each type can be used by most disci^^lines even 
though one cannot easily interchange lecture halls and offices. The assump- 
tion of continuously variable space is a weakness of the model, because new 
construction occurs by project or building and, therefore, occurs in quan- 
tum jumps. However, we do include the time lag of construction from start 
to completion. Recalling that ^12^^^ available space at t , we let 

k^ be a scalar weight and 

^l'^2'^3 space standard coefficients determining space needs 

as a linear combination of total faculty (TF ) and total students (TS ) ^ 
The fifth term of is then 



Student Mix 

Fiscal planning can be much more effective if student enrollments 
can be forecasted several years into the future. SCPM does not attempt 
to describe enrollments by discipline, although it could by defining 
additional state variables and equations; instead, SCPM focuses on student 
levels. Furthermore, for small colleges it was felt that average costs 
would not vary significantly across disciplines because these small colleges 
rarely can afford massive commitments of dollars to facilities and 




(iii) 
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faculties in the hard sciences and engineering, traditionally the most ex- 
pensive disciplines. (Schools such as Cal Tech and MIT are clearly 
exceptions to this .-^''e and they would need to recognize student discip- 
line and level.) The next three terms of are constructed from the 
same pattern as for faculty. 

^.U' Sc(t) - S.(t)} , i = 6,7,8 , and S. students of level i . 

(iv) 

Monetary Balance 



The last term of is a balance equation, expressed as a quadratic 

penalty function, which forces the annual net cash balance at the end of 
each period towards zero. Campus administrators are assumed to seek 
policies so that the cash inflow, e.g., transfers from endowment, gifts, 
and student revenues, is equal to the outflow, e.g., transfers to endowment, 
faculty salaries, construction costs, maintenance costs, and other operating 
costs. Otherwise, too much is withdrawn from income producing investments 
or, conversely, not enough is invested — both of which have an opportunity 
cost to the institution and should be avoided. The net cash balance is 
given by x^^(t'fl) ; however it is included in as a function of 

x(t), u(t), z(t) . With 

kg = a scalar weight and 

^1/ • > SiA n defined as before, the final expression in V is 

12 



. (v) 
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Terminal Conditions 

Because this optimal decision problem is formulated as a finite horizon 
differential dynamic programming problem, it is necessary to introduce "ar- 
tificial'' targets in the last planning period to correct for the truncated 
horizon. To prevent SCPM from "selling off" uncompleted space in the last 
three periods (space which is not usable during the model's lifetime, and 
therefore, space of no va^ e to the ins::itution) we drive Xg (N) ,x^q(N) ,x-j^-j^ (N) 
to zero by including in P(N) 

2 

k^QX^(N) , i = 9,10,11 , where k^Q is a scalar weight. (vi) 

The last term of P(N) is the net cash balance equation expressed as a func^ 
tion of x(N) , and write 

2 

kgXj^^(N) , kg a scalar weight. (vii) 

When put in the form of equation (7) , expressions (i) - (vii) constitute the 
criterion function currently being investigated. 

As discussed previously, the scalar weights were chosen to equalize unit 
losses and thereby avoid inducing artificial minima. Table 3 shows the values 
for the scalar weights used in this study. 

Although the target ratios were chosen to be constant over time, it 
would be possible to choose target paths showing the evolution of the var- 
iables over time and SCPM would then solve for the optimal decisions for this 
path. Another formulation of the criterion function would enable the campus 
planner to determine the shortest time in which the desired targets could be 
achieved. However, the numerical results reported in the next section em- 
ployed equations (i) through (vii). 
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SUMMARY OF RESULTS 



To assess the realism of the model, a reference set of input data was 
developed which represents the operations of an "average*' small liberal art 
college. Actually, these data were modelled after data in the University 
of Santa Clara [1970]. Changes in this reference set (denoted Data Set 1) 
were made to reflect different operating policy decisions which might be 
made by campus administrators and to reflect different assumptions about 
future levels of external financial support. The optimal decisions for 
these new environments were then compared to the base case results. (See 
Tables 2 and 3 for a summary of Data Set 1). 

Additionally, the college's optimal decision problem was solved for 
only non-negative values for the control variables. The resulting values 
thereby achieve sub-optimal solutions which trade optimally for more rea- 
listic results. It was shown that these solutions achieve the targeted 
ratios and provide adequate space as specified by the objective function 
but that they result in a much greater variation in the net cash balances 
at the end of each operating period. (See Table 5 for a comparison of net 
cash balances.) 

Student admissions, faculty hiring, and new construction decisions 
are the normal outputs of the model. In addition, the model can answer 
questions about the operational effects of alternative funding methods 
which are of considerable interest to external funding agencies such as 
state and local governments. To examine these questions, three additional 
data sets were analyzed from this policy perspective. 
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TABLE 2 



Accounting Summary for SCPM Example* 
(Enrollment Level : 800) 



EXPENSE 


% of 
Total 


INCOME 


% of 
Total 


Instruction 


$ 


104,400 


5.39 


Tuition & Fees 


$ 970,400 


58.07 


Student Aid 
Student Services 

Total Student Related 
Salaries: Teaching Faculty 
Salaries: Admin. /Other Fac. 


$ 


119,200 
266,800 
490,400 
418,000 
222,000 


7.31 
16.33 
30.03 
25.60 
13.60 


Transfers from 
Endowment 

Gifts & Grants 

Other 


115,200 

364,000 
221 ,600 


6.89 

21 .78 

13..26 


$ 


(gov't, aid) 
TOTAL 


$1 ,671 ,200 


100.00 


Total Fac. /Admin. Salaries 


640,000 


39.20 








Plant M&O 




208,000 


12.73 








Plant Additions 




294,450 


18.04 








Total Plant 


$ 


502,450 


30.77 








TOTAL 


$lj632,850 


100.00 


: 







★ 

This accounting summary was derived from "An Introduction to Program Planning 
Budgeting and Evaluation for Colleges and Universities" - July 1970 - University 
of Santa Clara— Office of Institutional Planning. 

Additional Assumptions: 

1) Eighty percent of instructional costs are faculty salaries; 

2) There are $3,840,000 of restricted funds yielding 3% per annum; 

3) Plant M&O costs $2/ASF; there are 130 ASF/student; 

4) Student/faculty ratio is 16/1; faculty incl'ides teaching, research and 
administrative staff exclusive of support. 
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TABLE 3 
Data Set 1 



F-MATRIX 



10 



11 



12 



13 



14 



0.93 0.05 




0.85 0.05 




0.80 0.50 




0.10 












0.95 0.02 






0.70 0.03 






0.70 0.04 





1 

2 
3 
4 
5 
5 
7 
8 
9 
10 
11 
12 
13 
14 



1 .00 



1 .00 



1.00 0.98 



1 .00 



■1.7 -1.4 -1.1 -0.8 0.06 0.06 0.06 0.06 -4.5 -4.5 -2.0 -0.2 0.03 "1.00 



G-MATRIX 
1 c 



ERIC 



1 

2 
3 
4 
5 
6 
7 
8 
9 
10 
11 
12 

13 
14 



0.70 




0.70 




0.95 0.30 




0.10 






0.03 






0.60 0.02 






0.70 0.03 






0.70 0.04 






1 .00 










-1.7 -1.4 -1.1 -0.8- 0.06 0.06 0.06 0.06 -1.5 
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TABLE 3 (continued) 

Vector of Initial x-V alues 

x(1) x(2) x(3) x(4) x(5) x(6) x(7) x(8) x(9) x(lG) x(11) x(12) x(13) xjU ) 

10. 12. 15. 1.0 10. 170. 140. 130. 2.08 2.08 2.08 104. 384. 0.0 



Vector of Predetermined Funds ($10,000's) 



Per. 


Amount 


Per. 


Amount 


Per. 


Amount 


Per. 


Amount 


1) 


58.56 


6) 


61 .55 


n) 


64.69 


16) 


67.99 


2) 


59.15 


7) 


62.17 


12) 


65.34 


17) 


68.67 


3) 


59.74 


8) 


62.79 


13) 


65.99 


18) 


69.36 


4) 


60.34 


9) 


63.42 


14) 


66.65 


19) 


70.05 


5) 


60.94 


10) 


64.05 


15) 


67.32 







TARGETS AND WEIGHTS 



Item 


Target 


Weight 


Item 


Target 


Weight 


Assoc/ Ful 1 : 


1.5' 


1 .0 


Soph/Frosh: 


0.7 


T.O 


Asst/Full: 


2.0 


1.0 


Fr ./Frosh: 


0.65 


1 .0 


Inst/Full : 


0.5 


1.0 


Sr./Frosh: 


0.5 


1 .0 


Stu/ Faculty: 


16.0 


0.01 


ASF/Student: 


130. 


50.0 . 


Net Cash Bal: 


0.0 


0.5 
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Comparison of Optimal versus Sub-Optimal Results 



The general 



form of the optimal control problem is to 



Min 



u(t)l 



'J 



N~l 

= P[x(N)} + I V(;x(t),u(0,2(t),t} 
t=0 



subject to 



x(0) : 
x(t+l) 



known 



= f (x(t),u(t),z(t),tj . 



This formulation does not constrain' the sign of either the state or con- 
trol variables. This is the form for which the solution algorithm given i 
the Appendix was designed p While firing full professors and selling 
newly constructed space may appeal to some interests, negative values for 
control variables are in general not meaningful. Furthermore, our study 
revealed that unconstrained solutions to the optimal policy problem usual! 
contained several such negative decision values, which, while small 
in magnitude, were nevertheless inappropriate and unrealistic. 

•Rather than impose inequality constraints which would require a re- 
formulation of the model or attach penalty functions to the criterion 
function to facilitate the use of a sequential optimization algorithm, we 
initially included a switch in the computer program which set negative 

values of the computed control variables to zero within the iteration 
13 

sequence. Since the algorithm computes the improved u(t) values 
based on small variations in x(t) and since those variations are not 
substantially altered by zeroing out negative values of u(t~l) , the 

D 

One can solve the general formulation vjith inequality constraints 
at considerable additional complexity and expense, see Jacobson [1969], 
but we ^concluded that at this stage of development these refinements were 
not worth the additional cost because the constrained and unconstrained 
results were so similar. 
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resulting sub-optimal path should provide acceptable solutions. A compar 
Son of the two sets of results confirms this assertion. 

Three sets of comparison runs were made to determine the loss from 
sub-optimization. These runs were based on the original Data Set 1, Data 
Set 2 in which exogenous funds grow at 2% per year, and Data Set 3 in 
which per-student income increased by $100 per year. Tables 4 and 5 shov 
the differences in control variables and in the yearly net cash balances 
for Data Set 1. Similar results obtained for Data Sets 2 and 3. The 
differences between the two sets of control variable solutions for Data 



Set 1 are summarized below. 









Mean of The 




Mean Value 


Mean Value 


Absolute 


Control Variable 


Unconstrained 


Constrained 


Differences 


Full Professors 


0.12 


0.15 


0.25 


Assoc. Professors 


1.86 


1.85 


0.63 


Asst. Professors 


2.26 


2.35 


0.73 


Instructors 


4.74 


4.88 


0.30 


Freshmen 


286.39 


287.52 


5.03 


Sophomores 


23.37 


23.75 


5.39 


Juniors 


42.22 


42.62 


4.73 


Seniors 


8.12 


8.36 


4.01 


New Construction 


2.92 


2.94 


3.41 



(000' s ASF) 

It is clear from an examination of Table 6 that in terms of the tar- 
geted values specified by the objective function, there is very little 
loss of utility associated with using the constrained formulation. 
Because the negati\i;.e values for the control variables generated by the 
unconstrained solutions are unrealistic, we shall concentrate on the 
constrained, sub-optimal results for the balance of this discussion. 
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TABLE 4 

Comparison of Control Variables 
Data Set 1 
Unconstrained vs Non-Negative Modes 





NEW FACULTY 


Unconstrained 


Non-Negative 


Per. 


Prof. 


Assoc. 


Asst. 


Inst. 


Prof. 


Assoc . 


Asst. 


Inst. 


1 


0 1 D 


3 20 


5 11 


3.82 


0 10 


3 .20 


5.11 


3.82 


2 


-0.03 


1 .68 


1.43 


4!63 


0.00 


1.68 


1.43 


4.63 


3 


0.12 


1 .75 


2.25 


4 .43 


0.00 


1 .48 


1 0 


4.66 


4 


0.02 


1 .54 


2.03 


4.53 


0.16 


1 .89 


d.l9 


4.37 


5 


0.39 


2.12 


2.78 


4.67 


0.31 


1 .94 


2.67 


4.77 


6 


-0.21 


1 .23 


1 .36 


4.47 


0.00 


1.14 


1 .50 


4.85 


7 


0.43 


2.20 


3.13 


4.73 


0.12 


2.08 


2.45 


4.38 


8 


-0.06 


1 .46 


1.62 


4.63 


0.32 


2.00 


2.78 


4.93 


9 


0.22 


1 .92 


2.60 


4.66 


0.00 


0.38 


1 .08 


5.68 


10 


0.10 


1 .71 


2.21 


4.80 


0.00 


1 .38 


1 .77 


5.24 


11 


0.18 


1 .90 


2.36 


4.69 


0.97 


4.01 


4.18 


4.49 


12 


0.05 


1 .66 


2.18 


4.79 


0.00 


1.14 


1 .33 


5.14 


13 


0.25 


1 .98 


2.62 


4.90 


0.00 


0.62 


1 .53 


5.50 


14 


0.12 


1 .84 


2.14 


4.72 


0.00 


1 .73 


2.17 


5.18 


15 


0.00 


1 .56 


2.28 


5.02 


0.53 


3.25 


3.64 


4.95 


16 


0.53 


2.53 


3.07 


4.89 


0.00 


1 .05 


1 .36 


5.47 


17 


-0.42 


0.96 


1.13 


4.83 


0.00 


1.70 


2.05 


4.88 


18 


0.77 


2.81 


4.12 


5.33 


0.39 


2.67 


3.40 


5.04 


19 


-0.27 


1.37 


0.87 


4.56 


0.02 


1 .83 


1 .89 


4.77 





NEW STUDENTS 


Unconstrained 


Non-Negative 


Per. 


Frosh 


Soph. 


Jr. 


Sr. 


Frosh 


Soph. 


Jr. 


Sr. 


1 


270.8 


26.5 


42.5 


10.4 


270.7 


26.5 


42.5 


10.4 


2 


272.8 


20.7 


39.5 


7.1 


272.7 


20.8 


39.6 


7.1 


3 


272.8 


21 .4 


39.5 


7.1 


272.8 


21 .4 


39.5 


7.1 


4 


272.8 


21.3 


39.5 


7.0 


272.8 


21 .3 


39.5 


7.0 


5 


282.1 


27.9 


45.6 


11.7 


282.4 


28.0 


45.7 


11.8 


6 1 


271 .9 


15.2 


34.4 


2.3 


276.3 


18.1 


37.1 


4.4 


7 


285.5 


30.5 


48.0 


13.7 


277.6 


22.3 


40.7 


7.7 


8 


279.4 


18.4 


37.6 


4.5 


289.1 


29.8 


47.7 


12.9 


9 


284.3 


25.1 


43.5 


9.5 


282.8 


18.5 


37.9 


4.5 


10 


286.3 


23.8 


42.5 


8.4 


278.9 


19.2 


38.3 


5.2 


11 


287.9 


23.5 


42.5 


8.2 


301 .4 


37.4 


54.9 


18.4 


12 


288.1 


22.7 


41 .8 


7.5 


294.5 


19.4 


39.5 


4.6 


13 


294.1 


26.7 


45 5 


10.4 


288.6 


18.7 


38.7 


4.6 


14 


292.2 


21 .8 


41 .4 


6.7 


289.1 


22.7 


42.0 


7.6 


15 


294.5 


24.4 


43.8 


8.7 


307.0 


35.2 


53.5 


16.4 


16 


304.2 


29.8 


■ 48.9 


12-. 5 


300.2 


19.7 


40.3 


4.7 


17 


289.5 


13.8 


34.6 


0.6 


294.2 


19.1 


39.5 


4.6 


18 


315.3 


40.2 


58.4 


20.1 


309.3 


33.4 


52.3 


15.1 


19 


294.0 


10.4 


32.2 


-2.2 


302.5 


19.8 


40.6 


4.8 



TABLE 5 

Comparison of Unrestricted Funds Balances 

Unconstrained vs Non-Negative Modes 
(Data Set 1 ) 



Per • 


Unconstrai ned 


Non-Negative 


1 


$ 15,240. 


$ 15,220. 


2 


-24,096. 


- 47,390. 


3 


13,195. 


16,260. 


4 


- 5,572. 


26,660. 


5 


856. 


54 , 500 . 


6 


787. 


14,030. 


7 


678. 


45,840. 


8 


- 4,560. 


-139, o^U. 


9 


9 , 1 66 . 


- 97,000. 


10 


-11 ,996. 


146,450. 


11 


10,543. 


109,970. 


12 


- 3,486. 


-128,660. 


13 


- 8,290. 


- 52,060. 


14 


20,927. 


101,7,0. 


15 


-28,550. 


23,120. 


16 


25,150. 


- 11,670. 


17 


- 7,584. 


16,080. 


18 


-21 ,526. 


7,077. 


19 


25,923. 


14,157. 


NET 


$ 69833. 


$ 5,974. 
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Comparison of Control Variables 

Unconstrained vs Non-negative Modes 
(Data Set 1) 

New Construction 




10 12 14 16 
Planning Periods 
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Comparison of Results from Various Policy Alternatives 

Perhaps the most striking feature of the solutions generated by the 
model is that they do not show a smooth expansion path, either in terms of' 
total enrollments or in terms of any of the control variables. This is not 
surprising mathematically, but it may surprise administrators unfamiliar wi 
controllable, dynamic systems with different response times. 

Intuitively, what has happened to the reference data set for the con- 
strained case is the following. The preference function desires the ratio 
of available ASF per student to be fixed at 130. Enrollments are therefore 
constrained during the initial four years by the amount of physical space 
under construction at the beginning of the planning period, (^g»^;^Q» x^ 
at t = 1) . For the reference data set, these are set at a level which ex- 
actly counteracts the depreciation of the existing capital stock. For the 
first four years therefore, the physical space available remains constant. 
Because faculty is linked to enrollments through the Student/Faculty Ratio, 
it is unnecessary to spend operating dollar balances resulting from the 1% 
growth in outs' funds, to increase the size of the faculty during this 
period. As a resuxt, the model spends any '^excess funds'* on new c-qpital 
construction. This becomes available after the fourth year at which time 
enrollments begin to increase. 

Although students yield a net dollar gain from the tuition level, they 
induce costs in the form of faculty salaries, capital needs, etc., so that 
it is impossible to pay for additional space by simpl , adding students to 
the rolls. When enrolliii::^nts increase, therefore, the amount of new con- 
struction must decline relatively. This will eventually cause some decline 
in enrollment levels from peak periods, and thus induce a cyclical pattern 
of physical expansion and student growth. 



41 



Although the system dynamics and the solution process are considerably 
more complex, the foregoing is the predominant reason for the cyclical var- 
iability in the expansion path. Another reason is the different time con- 
stants or response times for the various state variables. Full professors 
spend a longer average time in the system than instructors, freshmen more 
than entering seniors, and so forth. The time behavior of aggregates is 
built up from many of these overlapping transients and, therefore, the ag- 
gregates show a cyclical time behavior. 

Data Sets 1 and 2 differ only in the assumptions concerning the levels 
of external aid, with a 1% yearly increase in gifts, grants and government 
assistance reflected in Data Set 1 and a 2% yearly increase contained in 
Data Set 2. Both employ the same levels of aid for the initial year. At 
enrollment levels of 800 students, these funds account for 35% of current 
income in the initial year and 36.02% and 37.91% respectively for the total 
planning horizon. A comparison of the enrollment patterns generated is 
shown in Table 7. As can be seen from the table, there is no substantial 
increase in enrollments until the 10th period despite the fact that an 
additional $225,700 has been received by the institution through the 9th 
period. Since the 19-year "marginal-cost" of increased enrollments is ap- 
proximately $5,500, one might intuitively expect that a smooth expansion 
path would be generated which should have enrolled an additional 41 stu- 
dents through the first 10 years. In fact, the school has only enrolled 
an additional 3 students. 

Another way of looking at the data is to examine the increase in 
external funds in relation to the total income over the 19-year period. 
Here it can be seen that the additional $1.1876 million generated by Data 
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Set 2 is 3.51% of the total income generated via Data Set 1. This 3-5% 
increase in total income thus would result in a 1.35% increase in total 
enrollments if it were to come in the form of additional gifts, grants, 
and government aid according to this exponentially increasing function 
of time. 

Suppose instead the increase in income had come from tuition charges. 

Data Set 3 simulates the effect of a $100 increase in tuition charges (or 

equivalently , an additional subsidy of $100 per student) ^y changing the 

net return from students from $600 to $700. (Recall that the coefficients 

f ^ , . and / . , 1=5-8 , are the difference between tuition and average 
14, J 14, J 

student-related costs such as adxaissions, counseling, student-aid and health 
facilities . ) 

If the institution were to increase tuition and fees by $100, neglect- 
ing per-student student aid increases, enrollments would increase 3.57% 
over the 19-year period. The increase in total income represented by this 
policy change is $2,356,449 or approximately 7%. Table 8 shows a compari- 
son of enrollments and income under the assumptions of Data Sets 1 and 3. 

Looking at these results from the point of view of a potential fundor, 
such as the State or Federal government, it would appear at first blush 
that the most productive means of funding the inscitution would be the per- 
student instituticrial subsidy. In order to investigate the question more 
thoroughly. Data Set 4 was developed. This data set increases the external 
inputs by an amount equal to the $100/student-additional income generated 
by the increase in tuition at optimal enrollment levels. Table 9 shows 
the enrollments generated by Data Sets 2, 3 and 4 and the marginal costs 
of enrollments. A comparison of these costs shows that for the purpose of 
increasing enrollments, it is more effective to fund the institution directly 
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rather than paying a per-student amount. The same subsidy, in terms of 

both total dollars and timing, has been assumed by Data Sets 3 and 4, yet 

the solution to the planning problc. generates greater enrollments when 

14 

the subsidy is given as a fl^t grant. 
Conclusion 

These examples show that realistic and relevant results can be obtained 
relatively easily and inexpensively by SCPM. The control theoretic approach 
both incorporates the multi-level, multi-decision maker hierarchical struc- 
tures of higher education and enables educational planners to derive improved 
institutional plans and to evaluate many alternative operating policies. The 
robustness and flexibility of SCPM suggest that it could make a major con- 
tribution to improved educational planning. 



This counterintuitive result does not hold for the unconstrained 
case, if the optimal value of the criterion function is zero. If it is 
possible to achieve all the targets exactly in all planning periods, then 
MIN J = 0 and the optimal enrollments for Data Sets 3 and 4 would be 
identical. An alternative view of student tuition models and the effects 
of gov'rxnment subsidies is given in Weathersby [1970]. The above conclu- 
sion is supported in this supply and demand analysis. 
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APPENDIX 

Solution Procedure 

The linear-quadratic minimization problem is solved in one step using 
an adaptation of an algorithm devised by David Mayne [1966]. The following 
computer program is based on this algorithm. 

A non-optimal trajectory is generated using a nominal control sequence. 
The effect on the criterion (penalty) function of small variations of the 
control sequence is determined. This enables an improved sequence to be 
chosen. In the case of a quadratic criterion and linear system dynamics, 
the first improved sequence is optimal. 

The advantage of Mayne* s approach over conventional dynamic programming 
approaches lies in immense reduction of core requirements. In place of the 
optimal return function V° of Bellman, Mayne uses VV° , the optimal 
variation in the non-optimal return function due to variation of the state 
variable • VV° is expanded in a power series (to second order) of the 
variation in 'x , and difference equations are derived for the coefficients 
of the series. It is the identification of these coefficients which pro- 
vide, in analytic form, the optimal change in x(t) , and hence, by working 
backward in time, of u(t - 1). 

The program, in its present form handling 25 variables and 20 planning 
periods, requires approximately 50K bytes of core and 1,5 minutes CPU time 
on IBM's 360-65 O.S, On the University of California Administrative Data 
Processing System this costs approximately $10.00 for an adequate numerical 
solution for one case- It is hoped that future versions will reduce the 
^ize and cost of running this program- 

ERIC 
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COMPUTER PROGRAM 



FORTRAN IV G LEVEL 



18 



MAIN 



DATE « 71300 



l0/27/t>9 



0001 



0002 
0003 



0004 
0005 
0006 



THIS PKOORAy. EMPLOYS AN ITERAriVE PROCfcJURE TO FIND A SEQUCNCE 
3F CCNTROL VECTJHS WHlC-i MINIMIZES A NON-LiNEAR ^4-PER:0U OBJECTIVE 
FUNCTION SiGJECT TO LINEAR-DYNAMIC CONSTRAINTS. THE Pr<OGRAM IS 
3ASED 3N AW ALGORITHM DEVELOPED BY DAVID MAYNE, PUBLISHED BY 
INT. •J3URNi^L DF CONTROL < 1966). 



COMMON X(14,20),UC09,19) ,.MXf NUtNiPP .VX( l^f ) , 
1 VUJ(9, 9) ,VXX2( i'Vf 1^ ),VU( 9) ,LPP, VX2( l^) ,V( ?.0) f 



2: lN*iOUT ,NUMRUN 
RJ=AL JJ.OPT 



VXX(l^,l^),VUX(9.U>, 
F ( 1^, ,G( l<r,9) , 



RJ=AL JJ.OPT 

DIfJENS ION A( , B(.9, ) , AL ?HA ( 9 f 19 ) ,BETA( 9t I^, 19) U( 9) t 

lH(l<rfl4) .Z(l^tl9)fXNEW(X4f20),DELU(9)tW(l<rf l^)f 
2C(9,9),ClNV(9.9iffCC(A5l 

DCUSLE PRECISION C#CINV » CC 

I IN^l 

ICUT=6 



O0G7 
0003 
0009 
0010 
OCll 
0JI2 
0013 
001^ 
0015 
0016 
0017 
0013 
0019 
0020 
0021 
0022 



0023 
0024 
0025 
0026 
0027 
002 8 
0029 



00 30 
0C31 
0032 
0033 
0034 
0035 
0036 
0037 
0038 



^EAD BASIC DATA 

1 REAOI IIN,900) NLlMIT.LPRNT.NUfNX.f^PPf NUMRUN.OPT 
IF ( NLJMIT,GT.99) STOP 

LP? = NPP+1 

Ir(N'jMr^.UN.GT.l) GO TO 150 

2 DC LOO I^ltNK 

IDO READ (I lN»901)(Fa* J )t J^lfMXI 
I FCNUMRUN.GT.l ) GO TO 150 

3 00 101 I = l»SiX 

131 READ( I IN,90l ) (G( I, J ) t J^l.NU ) 
IF(WUMRUN.GT.l) GO TO 150 

4 DO 102 I = l,r»iX 

.102 REAO( I lNf90l )<H(I, J>, J«ltMX) 
IFINUMRUN.GT, I) GO TO 150 

5 DC 103 l«l.NX 

103 READ(IIN,901) (Z(l,J),J«ltNPP) 
IF(NUMRUN.GT»l> GO TO 150 



: INITIALIZE SYSTEM 

6 R£AD(UN»901)(xn.I)rI«ifNX) 
I F( NiUMRUN.GT .1 ) GO* TO 150 

7 DC 110 I»l,NU 

ILO REAO( I1N.901) (U(l. J) tJ=*ltNPP) 

IF( NUMRUN-EO. 1» GO TO 1000 
150 READ( I IN.922) ICHMGE 

CO TO (I000t2,3t4,5,6,7) t ICHMGE 

: CALCULATE INITIAL STATE VARIABLES 

1600 NCYCLE«0 

DC 120 L«1,NPP 
H^L-1 

DO 121 ;=1,NX 

121 X(ItM)=0.0 

DC 120 1=1 .NX 
00 122 J=l»NX 

122 X (I »M) «jX( I ,M)^F( I, J)*X(J»L)*H( If J ) ^'ZCJtD 
DO 120 J=1,NU 
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0039 



OC^O 
0041 
0042 
0043 
0044 
0045 
004S 
0047 
0043 
GQ49 
0050 
0051 
0052 
0053 
0D54 
0055 



00b6 
0057 
0C5S 
0C59 
0'./6C 

o;oi 

0062 
0063 
0064 
0065 



0C66 
0067 
0C68 
0C69 
0070 
0071 



0072 
0073 
0G74 
0075 
0076 
0077 
0078 



0079 
0060 
C08I 



120 X(I,M)«X(I ,M)+G(I,J)«U(JtLI 
ECHO CHGCK initial DATA 



WRirE( rQUT,902) MUMRUN 

W«Jrt< lOUTf 903> 

DO ?00 1=1 f NX 
200 WRITE( I0UT,904> 

WRITEC I0UT,905) 

WRITE( I0UT,903) 

DO 201 1=1, NX 
231 WRITE( I0UT,904) 

WR!Te(IGUT,906) 

WRire( I0UT,903 M 

DO 202 1=1 ,NX 

202 »JRITt5( iaLiT,904) 
WRITE( I0lJT,907 ) 
W?.ITE( IDUT,903) 
DC 203 J=1,NPP 

203 WRI TE ( I0UT,904) J 



, 1=1, NX) 
( F( I ,J) , J-1,NXI 
,I=>=1,MU) 
(G(I,J) ,J»l9NU) 
I«1,NX) 

(H( I ,J) ,J»1,NX) 
,i = l.»NX) 



C 
C 
C 



C 

c 



( Z( I ,J I , I»l, KX) 

BEGIN IVERATIVE PRDC^DURE 

CALL CALCV(O) 

300 WRITE( I0UT,908) NSYCLE 

WRITci ICUT,909) ( l»l=l,NJ) 

DC 301 J = l ,NPP 

331 WRITe( 10UT,910) J , ( U( I , J ) r I « 1 f hiU ) 

WRITE { [OUT ,911 ) 

VJRITEC IOUT,909) ( J, J«1,LPP) 

DO 302 r^l ,NX 

302 WRI Te(lOUT,912) I » ( Xi I » J ) » J*^l , LP P) 

CALL CALCV(LPP) 

CALCULATE J-VALUE 
J J=0.0 

00 310 I»=1,LPP 
310 JJ'^JJ+VI I ) 

WRITlK I0UT,913 M I , 1=1 .LPP) 
WRITE ( IQUT,914) ( V( 1) « 1 "1 ,LPP) 
lJftiT£( I0UT,915)J J 



:HECK TO SEE IF CURRENT CONTROLS ARE OPTICAL 

CHECK-ABSC JJ-OPT J 
IF(:HGCi<,Gr*0.02) CO TO 320 
WR;rE( I0UT.916) 
GO TO 1 

320 I F::^JUKIT--^4CYCLe) 321, 322,322 
3>1 V^^^ITiK lOUT ,917 ) 
GO TO 1 



C 

c 
c 



3EGIN MAJOR LOOP 

322 AA»0-0 

00 400 n«l,N?P 
I T=LPP-I I 
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0082 



0122 
0123 
Ol2<r 
0125 
0126 
0127 



C 
C 
C 



CALL CALCVUT ) 
CALCULATE MATRIX / 



VXX + F • VXX( T*^l ) F 



0033 






00 401 J=l,NX 


0084 






00 401 1=1, MX 


0085 






A{J,J)«:0.0 


00 36 






DO 401 K=l ,NX 


00^7 




401 


A(I,J>=A(I,J)+F(K, U*VXX2(K, J) 


00 






DC 402 J^lfNX 


00&9 






00 402 1=1, NX 


0090 






W II. J 1=0.0 


00<J1 






00 ^02 K=l ,MX 


0092 • 




4D2 


W ( I , J ) « W I I , J ) ^A( I • K ) «=F ( Kf J > 


009 3 






DC 403 J«lf\'\ 


0094 






00 4C3 I=«l ,:\JX 


0095 


C 


403 


All, J)«:W( I,J )+VXX( I ,J) 






:ALCLLAT£ MATRIX 3 «= VUX + G' VXX(T + 1) 


0006 


c 




OC ^04 1=1, MU 


0097 






00 404 J = l t'^X 


C0-?& 






S( I, J)«=0.0 


0099 






DO ^04 K=:l, NX 


0100 




404 


S (I , J)«=D( I ,J)<-G(K,n*VXX2(Kf J) 


0101 






00 ^05 1=1, MU 


0102 






DC 405 J^-1 tN'X 


01U3 






r.d , J) ^0.0 


0104 






on 405 K«1,NX 


GIGS 




4:) 5 


ml , J>»W(I ,J)^6( I«K)«F(K, J) 


01^6 






DO 406 I«=1,MU 


0107 






00 ^06 J«1,NX 


0103 




40 6 


B (I ,J)<=W( 1 ,J)«-VUX( I ,J) 




c 
c 

r 


C^ICULATE MATRIX C » VUU + VXX(r+lJ 


0109 






DO ^07 I«=1,NU 


Olio 






DC 407 J=^1,^JX 


OUl 






W(It J>*'--0,0 


Olli 






DC 407 K=1,NX 


i) 1 1 '5 




40 7 


W (I,.J J^V.'( I, J )+G(K» i )<-VKX2(K,J) 


0114 






DO y08 I =1 ,NU 


0115 






00 403 J^lt'^U 


0116 






C{ I, J J =0.00 


0117 






00 40b K=l ,NX 


OllB 




408 


ClU JJsCl l,J )+W( I»K)*G(K, J) 


0U9 






DG 409 2=1 ,NU 


0120 






DO 409 J«l ,NU 


0121 




409 


C(I,JI»C<I,J)+VIJU<I,J) 



Z PLACE UPPER TRIANGLE OF C INTO CC BY COLS, 
C 

00 419 J-^l.NU 
OC 41B I«=1,J 
IF(C( :,Jn 415,416,416 
415 CC(< i^DMAXK C( I9 J ) tC( J,I ) ) 
GC TO 418 
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0126 ^16 CC(.<I«0MIN1<C( ItJ)tC(J, n> 

0129 418 K=K+1 

0X30 419 CONTINUE 

0131 IER = 0 

0132 PPS«*000001 

0133 CALL OSlNVCCNUtEPSiieR) 

■c 

C PLACE CC IN CINV 
C 

01.34 L-1 

0:l35 DC ^85 J»ltNU 

0136 00 480 I*ltJ 

0137 C INV( I , JJ^CCtL) 

0133 IFd.NE.J) CINVIJ,! )«CCa) 

0139 t^BO L-=L+1 

0140 48$ CCNTTKUE 

0141 IF( lER.EO.-ll GO TO 1 

0142 DO 410 J=1*NU 

0143 HUU>=0,0 

0144 DO 411 I-^l.-MX 

C145 411 H'J(j}=HUU)<-VX2(n*G(I,J» 

0146 410 HC!( J) =HU(J)^-VU( J) 

c 

C CALCULATE ALPHAdT > • B£TA( IT ) 
C 

0147 DO 420 I=»1.NU 

0148 ALPHA( I, IT )^0,0 
01^9 DO 421 J=l iNU 

0150 421 AL?HA< I, lT)::iALPHA( I , I T ) ^-C I MV( I »JJ*HU(a) 

GI51 420 ALPHAi I, IT 3=-ALPHA( I, ITJ 

0152 DO 430 1=1 tNU 

0153 DO 430 J=ltNX 

0154 BETA! I , J,lT)=0.0 

0155 . DO 431 K=ltMU 

0156 431 B£TAU,.J, IT )»B£r A( ItJ, ITKC INV( I »*BCKf J) 

0157 430 B€TA( I ,J,lTl=-eETA<I.J,IT} 
C 

C CALCULATE NEW VXX2 AND NEW VX2 
C 

0158 DO 440 I=1»MX 

0159 W(Itl)«0,CJ 

0160 DO 440 J-1 ffiU 

0161 440 W (1, 1 )-W( i;v 1) i-BETAC J, I tiT) *HU{ J) 

0162 DO 441 I = :,,NX 

0163 V^(I»2)=0.0 

0164 DO 441 JrltMX 

0165 441 W(I .2)«=W(I t2)>VX;M Ji*F(Jtn 

0166 . DO 442 I«1«NX 

0167 V42 VX2( n=VX< I)+W'iKtll*W(i,2) 

0168 DO I=ltMU 

0169 . 00 450 J=i,MX 

0170 W(I»JI=0,0 

0171 DO 450 K==l.NU 

0172 450 W (I.J )«W( I.J )^CINV( I,K)*B(KtJ) 

0173 DO 451 I=^ltNX 

0174 DO 451 J»l oNX 

0175 VXX2( ItJ)«0.0 

0176 DC 452 K«ltNU 
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0177 
0178 



0179 
OlSO 
*18I 
0182 



452 VXX2(I tJj^VXX?. (IiJ)+B(K,I)*W(K, J) 
451 VXX2(i[ tJ)=A< I, JJ-VXX2(I ,J) 



C CALCULATE AA TO INDICATE EXPECTED IMPROVEMENT IN J-VALUE 



W (If 1 ^''0 .0 
OC 460 I'^l ,NU 

460 w(i,i)«wn,n*Hi;(n*ALPHA(i, 

4D0 AA-AA+0.5*W(lf I) 



IT) 



:ALCL/LATE new controls and new STATE VA?=tIABt£S 



0163 
0 184 
01^5 
0186 
0187 
0183 
Oia<3 
0190 
G191 
0192 
0193 
0194 
0195 
01v6 
0197 
0193 
0199 
0200 
0201 
02 02 
0203 
0204 
0205 
0206 
0207 
0208 
0209 
0210 
0211 
0212 
0213 
0214 
021 5 
0216 
0217 
0218 
0219 
0220 
0221 
0222 
0223 

0224 
0225 
0226 
0227 



OC 500 I«l ,MJ 

U(l f 1>=U( I ,l)+ALPriA(I ,1) 

IFUd.D.LT.O.) Udf l) = 0. 
330 CONTINUE 

DO 501 I«1,NX 

XNEW( 1,2 )=0.0 

DO 502 J=l ,NX • 
50 2 XNlW( 1,2 ) = XNE^( 1,2 )+F(I . J) *X( Jf 1 J+H( I, J)*Z< J,l) 

DC 501 J=l ,NU 
501 XNEWd ,2)=XKRW(I,2)+G(I,JJ*U(Jtl) 

DO 510 L=2,NKP 

M = L*-1 

DO 511 I=>1 ,NU 

0£LU( I )=0.0 ■ 
DC 512 J«1,NX 

512 DELU(l) BETA( 1 ,J,L)«( ?Nr:W( J,L)-X( J.C) )*'0ELU(1 ) 
U(I.L)= OELUt I)+ALPHA(I,L)+jn,U) 

lF(Un tD.LT.a.) U(I»L) = 0. 
511 CONTINUE 

00 513 1=1, NX 

XNEWd ,M)^0.0 
■ 00 514 J=^1,NIX 
51V XNEW( I .M)=XNEW(I,M)+F(I, J)«XNEW( J,L>+H<I,J)«Z( J,L) 

00 513 J«l ,MU 

513 XS£W'{ I,M)«XNEW(I,M)+G( If J)*U(Jf L> 
510 CONTINUE 

DO 520 1=1. NX 

OC 320 J-?,LPP 
520 X(I , J) =XNEW( I f J) 

N CYCLE^NCYCLE*-! 

Wi^lTE{A0ljT»918)AA 

GO TO 300 
9D0 FORl-iATtoU.FS.O) 

901 rCTMAT (lOFa.O) 

902 FOf;MAT( 1H1.20X,«DATA FOR RUN NUMBER' • I3//20X, • F-MATRI X • ) 

903 FORMAT (/3Xrl4l9) 

904 FCRMAr(lX,r2.l4F9.2) ^ 

905 F0R^ATC//20X. •C-M^TR1X• > 
90 6 FCRV.AXxlHL ,20Xt • H-MATRIX • ) 

907 FCRMAiT(//20X.»Z-TRANSPOSE* ) 

908 FORMAT(lHi,lOX,«STATISTlCS FOR CYCLE NUMBER* t I 4//2 OX CONTROL VARl 
lAELESiM 

909 FijRMAt(/(3X»I0I12) ) 

910 FORMAT ax, 12, 5X, L0G12.4> 

911 FCRMATIIH1,20X,«STATE VARIABLES') 

912 FORMAT (/2X, 'Xt S 12, M • , 10C12. 4/ ( 8X, 10G12.4/) ) 
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DATE » 71300 



10/27/59 



PURPOSE 

INVERT A GIVEN SYMMETRIC POSITIVE DEFINITE MATRIX 



USAGE 

CALL DSINVI A.NfEPS, I£R) 



C 
C 
C 
C 

C 



DESCRIPTION OF P:;RAMETeRS 
A DOUBLE PRECISION UPPER TRIANGULAR PART OF GIVEV4 SYMMETRIC 
POSITIVE DEFINITE N BY N M^VTRIX, 

ON RETURN A CONTAINS ThE RESULTANT UPPER TRI ANGULA R MATRIX 
IN DOUBLE PRECISION, 
N ORDER OF THE GIVEN MATRIX 

EPS- SINGLE PRECISION INPUT P^RA^ETER WHICH IS USED AS RELATIVE 

TOLERANCE ri->K TEST ON LOSS OF S IGNI F ICANC E • 
I ER RESULT IMG F*^ROR PARAMETER CODED AS FCLLOh'S 



lER' 
lERx 



I cR = 5 



NO ERROR 

NO RESULT BECAUSE OF WRONG INPUT PARAMETER N OR 
BECAUSE SOME RAOICAIvlO IS NON-POSITIVE (MATRIX A IS 
IS NOT POSITIVE DEFINITE, POSSIBLY DUE TO LOSS OF 
SIGNI FICANCe, ) 

WARNING VJHiCH INDICATES LOSS OF SIGNIFICANCE. 



C REMARKS 

C ThE L'PPER TRIANGULAR PART OF GIVEN MATRIX IS ASSU^lEO TO BE STORED 

: CCLUMilWlSb IN rvj''(N*-l) /2 SUCCESSIVE STORAGE LOCATIONS. IN THE 

C SAME STORAGE LGCArlONS THE RESULTING UPPER TRIANGULAR MATRIX IS 

C STORED COLUMNWISE TOO. 

C TH6 PROCEDURE GIVES RcSULTS If N IS GREATER THAN 0 AND ALL 

C CALCULATED RADICANDS ARE POSITIVE. 



0001 

0002 
0003 



000^ 
0005 



0C06 
00C7 



CCOB 
0C09 
OOlO 
OOil 
0012 
0013 
0014 



C 
C 



C 

c 



SUBROUTINE REQUIRED - OMFSD 

SUaf^OUT INE DS!r>iV{A,N,ePS,IER> 

DI.SENSIGN Ad) 

DOUBLE PRECISION A, DIN, WORK 

FACTCRIZE GIVEN MATRIX BY MEANS OF SUBROUTINE DMFSO 

A=TRANSPOSE(T) ♦ T 
CALL UMFSDC A,Ni,EPS, leR) 
IFtlER) 9tl#l 

INVERT UPPER TRIANGULAR MATRIX T 
PRfiPARE INVERSION-LOOP 
I iPlV=\*{ N*l)/2 
IND==1PIV 

INITIALIZE I^4VERSIaN-L00P 
00 6 I»l ,N 
DIN=^1.DG/A{ IPIV) 
A( IP :v )«=DIN 
M!N-N 
KEND=I~1 
LANF«N~KENO 
IF(KEND)5,5,2 
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DSINV 



DATC 



71300 



lC/27/t>9 



0015 



0016 
0017 
0018 
0019 
0C20 



2 J-IND 

JNITIALIZE ROW-LOOP 
00 A K«l tKEND 
WORK=0.00 
M IN=M1N-1 
LHOR=£PI V 
LVeR=J 



0021 
0022 
0023 
0024 



START INNER LOOP 
DO 3 L = LANF.v^IM 
LVER=LVER+X 
LH0R = LH0R4 L 

WOKK=WDRK+ArLVERi*A(LHOR) 
ENO OL INNER LOOP 



0025 
0026 



A{J)°-W0RK*01N 
4 J»J-MiN 

END OF ROW-LOOP 



0027 
002b 



0029 
0G30 
0031 



0032 
0033 
0034 



0035 
0036 
CJ037 
003a 



0C3S- 
0040 



0041 
0042 



I P1V=IPIV- 
I N0«IND-1 
END OF 



MIN 

INVEP.SION-LDOP 



CALCULATE iNVERSEtA) 8Y MEANS OP INVERSElT). 

INVERSE! A i.--:INVeRSf ( T » *T RANS POS E ( INV ERS t ( T } ) 

INITIALIZE MULTIPHCATiON-tOOP 
DO S 1=1, N 
IPIV^^IPXV*-! 
J=«IPIV 

INITIALIZE ROW-LOOP 
DO 8 K^I .N 
•WCRX-O.DO 
LHOR«J 

START INNER LOOP 
DO 7 ..=K,N 
L VER«LKGR<-K-I 
WQR<=i<3ftK*'A(LH0RJ<-A(LVERI 
7 LHOR = LHOR'»-L 

END OF INNER LOOP 

A (JJ^woRK 

END OR ROW- AND MULTIPLICATION-LOOP 

9 RETURN 
ENO 
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DATE 



71300 



0001 



0002 
0003 



000^ 
0005 



0C06 
0C07 
0008 



C 
C 

c 



0Kf*m^t «<-'<i^)(i«ii4(4r«(4>««(4i ^ ^ i)i ifi ^i^^ ^ ^ 0^ 4:^1 ifL^f Kf^ ^ Kt Kt^t ¥ Kt ^ 9i 4^ <¥ Jtf tt^ ^U* ^ 

PURPOSE 

FACTOR A GIVEN SYM,^i^::TRIC POSITIVE DEFINITE MATRIX 
USAGK 

CALL 0MFS0( A,NtEPSt lER) 

DESCRIPTION OF PARAMETERS 
A U'JUBLE PK£;CISI0N UPPER TRIANGULAR PART CP GIVEN 

SYMMETRIC POSZTIVE DEFINITE N BY N MATRIK 
N THt NUKIUR OF POWS { COL UI^MS ) IN GIVEN MATRIX 

DN RETURN A CONTAINS THE RESULTANT UPPER 

TRIANGULAR fUTRIX IN DOUBLE PRcCISION 
EPS- SINGLE PHE^ISION X' n CONSTANT WHICH IS USED 

AS RELATIVE TOLERAwuE FOR TEST ON LOSS OF 

SIGNI FICANCE 
lER RESULTING ERRORPARAMETER CODED AS FOLLOWS 

IER»0 - NO- ERROR ' 

IER«-l - NO RESULT BECAUSE OF WRONG INPUT PARAMETER 
N OR BtCAUSE S0«1E RADICAND IS NON-POSITIVE 
(MATRIX A IS NOT POS IT I V EDEF IN I T E - 
POSSIBLY UUE TO L3SS OF SIGNIFICANCE) 

IER = K - WA«\JI^^G WHICH l^OiCMeS LOSS OF SIGNIFICAMCe 
THE RADICANO FORMED AT F ACTOR I ZAT I D\J 
STh'P K+1 WAS STILL POSITIVE BUT NO LONGCR 
GREATER THAN AS5 { E PS <'K -H ,:<+ Z ) ) 

REMARKS 

ThE UPPER. TRIANGULAR PART OF- GIVEN MATRIX IS ASSUMED TD tiE 
STDi^EO COLUMNWISE IN N«iNh-IJ/2 SUCCESSIVE STORAGE LOCATIONS. 
IN THE SAME SVOi^AGE LOCATIONS '.HE RESULTING UPPER TRIANGULAR 
MATRIX IS STORED COLUMNWISE TOO. 

THE PROCEDURE GIVES RC^SULTS IF N IS GREATLR THAN 0 AND ALL 
CALCULATED RAOICANDS ARE f^OSITIVE. 

THE PKOOUCT OF RETUR.'-^ED DIAGONAL TERMS IS EQUAL TO THE 
SQUARE-ROOT OF THE DrT£RMINANT OF THE GIVEN MATRIX. 

ME THOD 

SOLUTION IS DONE USING THE SQJARE-ROOT METHOD OF CHOLFSKY. 
THE GIVEN MATRIX IS REPRESENTED AS PRODUCT Of TWO TRIANGULAR 
MATRICES. THE LE, T HAND FACTOR IS THE TRANSPOSE OF THE 
THE RETURNED RIGHT HAND FACTOR. 

SUBROUTINE OHrSDI A,N,EPS,IEK> 



D IMcNSlON A( I) 

DCUSLE PRECISION DPIV,OSUMtA 

TEST ON WRONG INPUT PARAMETER N 
1 F(M-l )12»l ,1 
1 IcR = 0 

INITIALIZE DI AGNUNAL-LDOP 
KPI V=0 
00 11 K«!!.«N 
KFIV«KP)[ V^-K 
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DMFSD 



OAT£ s 71300 



10/27/59 



0009 
0010 



0011 



0C12 
0013 
0014 



0015 
0016 
0017 
001 B 



0019 
0020 



00 >l 
00>2 



•002 5 

0 0 ; * 
00 i 7 
00^8 



0029 
0030 



0031 
0032 
0033 
0034 



C 
C 



LeNO=K~l 

CALCULATE TOLERANCE 
. TCL=AB5(EPS*SNGL(A<KPIV) I ) 

START ^'ACTORI 2ATIQN-L0UP OVER K-TH ROW 
DO n I^K,N 
DSUX=0. DO 
i l=(LeND^/2t4i2 

START INMER LOOP. 

2 00 3 L=l .LEND 
LaNF = KPI V-L 

L ^^^o=^^Jo-L 

3 03U>;=0SjM + AtLANF >*A(LlNJj 

ENO Or INNER LOOP 

TFANSFORM ELEMCNT ' A C INO ) 

4 0 SUM*A( TNO) 

IFt l-K) 10»5»10 

* TEST FOR NEGATIVE PIVQT KLEHCNr AND FCK LOSS OF S IGNi I F ICAKC E 

5 Ir(SNGL(OSU^)-TDL> 6t6»9 

6 IFCOSUM) 12»12,7 

7 I Fl IFJt) B,3,9 
B IER=^K-l 

COMPUTE PIVOT ELEMENT 
9 OPIV^O^:OKT(DSUM) 
ACKPIVJ^DPIV 
C!?IV=1 .00/DPlV 

GO TO n 

CALCULATE TERKS IN ROW 

10 A ( INCJ~OSoM*DPIV 

11 iriO«lNDi-I 

LMO OF 0IA30t^AL LOOP 

RETURN 

12 !£R=-1 
RETURN 
6 NO 



ERIC 
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CALCV 



DATE 



71300 



10/27/59 



0001 



0002 



0003 
000^ 
COOS 
0006 
0007 



0008 
C009 
COlO 
COl I 
C012 



CO 13 
CO 14 
001*5 
0016 
0017 
0013 
0019 
0020 
0021 
0022 
0023 
C024 
002ti 
C026 
C027 
0028 
CC29 
C030 
0031 
0052 
0033 
0034 
003i> 
0055 
0037 
0038 
0039 
C040 
0041 
0042 
0043 



C 
C 

c 
c 



Suf--iQUT:.N5 CALCVdPFiR) 
IN WHICH IS CALCULAfC-D VX , V J • VXX i VUX • VUU, A\0 V 
THCfi£ ARt THP.tL SeCTlOKS TO THE ROJT INE 

1 PARAMETERS ARt READ IN, SECOND DERIVATIVE MATRICES SET TO 0,0 

2 V IS CALCULATED USlMG UPDATED X AND UiALSQ VX2(LPP) AND UXX2(LPP) 
SUM^iARY STATISTICS ARE PRINTED CUT 

3 VX,\/JfVXXi VUX, VUUf ARE CALCULATED FOR PERIOD IT AS DEFINED IN MAIN 

COMMON X( 14,20) ,U< 09119) ,NX, NU , NPP , VX ( 14 ) , VXX ( 14 . 14 ) • VUX ( 9 . 14 )# 
lVUJ(9t9 ) ,VXX2( I4i 14)» VU( 9) .LPPf VX21 14) , V( 20 ) f F ( 1 4. 14 ) . G( 1 4 , 9 ) , 

21 IN, lOUTfNUMRUN 
DIV.ENSIDN R(0),P(8) 
REAL K (9) 

iFdPER.EQ.O) GO TQ 100 
IF( IPE^.EQ.LPP > GO To 200 

GO TO 300 

READ BASIC DATA 



100 



103 



IF(NUMRUN.EQ. 1) CO TO 103 
READl I IN, 904) ICHNiGE 
T F(tCHNGE*EO.O) GO TO 101 
READJ X IN ,900) KiRiClvC2iC3 
WRITE! lOUT, 901 ) K,RvCl,C2.C3 



CALCULATE VXX.VUXiVUU (WHICH IS INVARIANT) 



DC 102 J«=l»14 
DO 102 1 = 1, J 
102 VXX(I.J)=0.0 

D=K( 1) *R( D-^t^a +K(5)'>C3**2 

VX>M 1)=:0+K(2 )'t'R(2)«*2+K( 3 ) *^ ( 3 ) **2+K ( 4 J ♦R ( 4) **2 

VXXC I,2)»0-K«2)*^>^(2) 

VXX( U3)«=D-K(3)*R(3> 

VXX(1,4)=:D-K(4)*R(4) 

V)(Xi2,2)=D + K(2 ) 

VXX(2,3) = D 

VXX(2,4)=0 

VXX(3,3) =D + K{3) 

VXX(3,4)«D 

VXXi4,4)=D+K(4> 

00 110 1 = 1.4 

DO 110 J=5i 0 

110 VXX( r » J) =~K(1 ) <«R<1 >+K<5)*C2*C3 
00 111 I»=l»4 

111 VXX( Ul2) = ~K{5)*C3 
D^'ti^i i ) +K('3 )*C2**2 

VXX( 5»» 5) '^O + KI 6)*R( 6)«'*2+K( 7) *i<< 7 )**2 + K(a )<*R( 8 ) **2 

VXX(5,6)=0-K(6)*R(6) 

VXX(5t7)?=D~K(7)»R(7) 

VXX( 5, a)=D~K{ 8)*R(8) 

VXX(6,6<=*0 + K(6) 

VXX(6,7)=D 

VXX (6, 8) ==0 

VXX(7r7)«D+K(7) 

VXXJ 7,8)=0 

VXX(8,8)«D*-K(8. ) 

DO 112 I«5f8 
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CALCV 



OATt = 71300 
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0C45 
0 046 
0C.V7 
004ei 
0049 
0050 

005: 

OC*jZ 
0053 
0054 
0055 
0056 
0057 

0059 
0060 
0061 
0062 
0063 
006^ 
0065 
0066 



0067 

OO'j^J 

GV/0 
0071 
CC72 
0073 
001<f 
0075 
0C76 
0077 
C07-J 
0 079 
OCtiO 
0031 

ooaz 

0063 
003^ 
00E5 
0086 
0067 
QOhH 
0069 
0090 
0091 
009.i 
0095 

009*V 
0095 

0096 



112 



121 



120 



130 



1? 1 



132 



13 3 



VXK( I ,12)« 
VXX( 12,12) 
DC 121 J«l 
OC 121 1=1 



-K<5)*C2 
= K{ 5J 
rl4 
» J 



= r;<9) fF( 14.n»F<l<».J) 



C 

c 
c 



VXX( I ,J) ! 
DC 120 
1 K1=I-1 

00 120 J^L, IMl 
VXXCi, J)«VXX(J,1) 
00 1 30 I «1 r8 
DG 130 j=i,a 
vuu( I «>n =vxxii , J ) 

VUX( I,J)=VXX(l.J) 
DC 131 [=l,tt 
DO 131 J=9,14 
VUX( I, J>«VXX( I, J) 
DC 132 J^^Ul-V 

VIjX( 9t J) «K( 9) *F(IV 9 J)*G{ 1** t9) 

DO 133 J=l,e 

VLUC9f J)=VUX(9,J) 

VUUt J,9I =VUU{9, Jl 

VCi;( 9>=K< 9)*G{ 14 1 9) «x<2 

GC TO SOO 



CALCULftTt PREFERENCe FUNCTlQNi AND VXX21 LPP) , VX2 ( LP!*) 



200 



202 



WRITE ( I GUT, 902) 
DO 210 J=L,N?P 
CG 201 I»l ,8 
P ( n -X( I , J)+U( 1 ,J) 
T 5=^0.0 

r F=o.o 

DO 202 1=1,4 
TF>-TF+P( I ) 
rS:=TSi-P( 1 + 4) 
P(2f=^>i2)/P(l> 
P£3)«=FM3>/P<1) 
P (^) =P(4) /P(l) 
P tl ) = T^:/TF 
A1 = ?16 >/P(5 ) 
A2=?< 7) /P( 5) 

A3=P<6)/P{5i ' 
A4=Ci ^-C2*TS+C3*TF 

WRITE ( IDUT,903) J,TStTF ,(P< 1) , I «=1 , 4 ) , Al , A2 t A3 , A4 , X (1 2 t J 1 

Al=RCl)*TF-TS 

A 2 = ?^ 'i 2 ) * ( X ( 1 , J ) +U t 1 r J n - C X ( 2 , J ) *-U ( 2 f J J ) 
A3 = ^{3i*(X{ 1,J) + U( l,jn-(Xt3,JJvj(3,J) ) 
A^ = R(£r)*lX ( I, J )+U( 1,J J )-(X(4, J >+U< 4. J n 
A5-=CUC2*TS+C3*TF-X{12, J) 

A6^?^ (6 >*tX< 5, J )Vj(5,J ) )-(X<6, J )+-U{ 6, JU 
A7 = .^(7)*(X<5,JUU(5,J) )-{X(7,J|+U(7, J)J 
A6=^^ ( a)*< X( 5» J )'»-U< 5,0 ) XI 8« J) -J-UC (i, jn 

V (J)=K ( L J«Al*>s^2+K { 2 )^A2**2<fC( 31*A3*v^ >K( 41 <«A4*« 2fK( ^J* A5«*2 
1+K(6)«A6*«=2 + K(7J*A7**2 + K(8 > * A3 **2+K{ 9 I *X ( 14, J + 1) *«2 
210 V(Jl=V{J)/2.0 

V<LPP)=0.5*K(9 )*X( l4,LPP)>i^*2 
1 ♦ 0. S"**;? ( 5) XC9,LPP)**2 + X(10,L^'f») ♦'*2 + X( 11 ,LPP)«»»2 ) 

DG 220 1-1,13 



ERLC 



I 
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0097 220 V>J?.{ I )=0.0 
0093 DG 2?A 1 = 1 , 1 <f 

0099 DO 221 J=l,l<r 

0100 221 VXX2( i , J »=«0.0 

0101 DO 222 1=9 ,11 

0102 V><2( 1 ) 5) fXi l,LPP> 

0103 222 VXX2{ [, I i^R( 5) 

0104 VX2(14 K(9) *X( 14,LPP) 

0105 VXX2( 14, 14)»=f^( 0) 

0106 GG TO 800 



C CALCULATE VX( fPER J 
C 

300 TF*=0.0 

DO .301 1«1,4 

301 T F = TF + X( I, IPER)+U( I, IPER I 

10 ^ X(5,rPER) + X«6,IPeR) +U(5f !PERJ+U«6, IPER) 
UO « X(7,IPeR) + XCe.IPER) +U( 7tlPER)+U< 8ilPER) 
TS » ..D + UD 
DO 302 I«l,14 
'302 VX( I J«0.0 
J=i PER 

A 1*-R( 1) t TF-TS 

A2 = R(25*(X( 1, J 5<-U(l,J ) )-(X(2f J )+U{ 2, J) I 
A2=R(3) 6^ (X(l, Ji + U< l,J) )-(X(3tJi^J<3, J) ) 
A-ii=:i(4)':^(X(l,J) + U<l,J)i-<X(4,J)+U(4,J)) 
A5=C1+C2«TS+C3*TF-X( 12, J) 

Ae=Rt6)«( X(5tJi+U(5,J) )-(X(o,J)*U(6,J) ) 
A7=;<(7)*tX(5,J)+U(5,JJ)-(X(7,J)+U(7,j)) 
A 6=:U 8) { X( 5, J) *-Ut^i, J) )-<X( e,J ) *-U(8, J) ) 

D="i\( 1)*A1*RI 1 )-»-K(5) <*A5*C3 

VXtl >=0'i-K(2 )*A2*R( 2 ) ( 3 ) * A3*R ( 3 ) + K ( 4) *A4«R ( 4) 
VX(2)=:D-K(2)v^2 
. VX(3) = 0-K( 3)*A3 
VX(Cr)=0-K(4)*A4 
D=>-K( l)»AH-K( 5J*A5*C2 

VX(5 ) = DfK( 6) <'A6*R( 6)+K(7)*A7*R( 7)+K( 8)*A8*R( 8) 

VX(7)^0-K( 7) 4'A7 
VX(8)=*0~K (8)*A8 
VX(12J=:-K(5)*A5 
00 303 1 = 1 . 14 

5D3 vxn )=vx ( n+K (9)*x( 14, j+i)*F( 14, n 

DO 310 1=1 .6 
310 VUi I )=.-VX( I ) 

ViJ(9)=.v(9)*X(l^ , J+l J«G( l4t9i 

800 COMTli\UE 
RET'JRN 

C 



0142 900 FQRKAT(9F8.0/eF8.0/3F8,0) 

0143 901 ^^CRKAr(lHl,20X«PREFERENCE FUNCTION WEIGHTS •// 1 X« KI-K9' , 9F 12, 6// 

1 1 X« kl-a8« ,8F12,6//1X,«C1~C3«,3F12.6) 

0144 932 F QRMAT (IHl , 50X, • SUKKAaV ST A Tl ST I CS • / / 1 X , • PE R • , 4X , "TS • , 9X , • TF • , 6X , 

1 ' TS/TF " , 6X, • F2/rl« , 6X, • F3/ K 1 • , 6K , • F4/F 1 • , 6X , « SG/F R • ,6X,« JR/FR' . 
26X,« SR/FR* ,5X.«ASF REQ» .6Xt«ASF BUILT*/) 

0145 9D3 FORMATC 1X,I3,2(3X,F 8. l),7(3X,r 8.3) •2(3X,F9t2)) 

0146 934 F CRViAT(ll) 

0147 END 



0107 
OlOa 
0109 
0110 
0111 
0112 
0113 
0114 

on 5 

0116 
0117 
C113 
0119 
0120 
0121 
0122 
C-12 3 
0124 
0125 
0126 
0127 
0128 
0129 
0130 
0131 
0132 
0133 
0134 
D135 
0136 
0137 
0138 
0139 

ri4o 

0141 
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DATA INPUTS FOR MULTIPLE RUNS 



To make multiple comparison runs without respecifying all the data 
inputs, coded data-change cards are used. Changes may be made in either 
the. data inputs to the main routine or to CALCV or both. 

Changes in the data for MAIN must be preceded by a RUN card and 
followed by a' card with "1" in column 1. Changes made in CALCV data must 
be preceded by a card with **2" in column 1. If no changes are made in 
MAIN data, a card with **1" in column 1 must follow the RUN card. If no 
changes are made in CALCV data, a card with "0*' (zero) must follow the 
last card in the data set which changes MAIN, i.e., the card with "1" in 
column 1. The END card which follows the data in single runs is to be 
made the last card of the change-data sets. (There is only one END card 
for the program.). 
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FORM OF INPUT DATA 



I. DATA FOR MAIN ROUTINE 



Item 


Var . Name 


Va 1 ue 


if LO 1 S . 


if KOWS 


rorma z 


1) RUN card: 










(l3,3x,4I3,F8.0) 


a) # Iterations this run 

b) # Control Variables 

c) # s^.tate Variables 

d) # Planning Periods 
c) Tvun Number 

f) Est. Optimal Value of p.f. 


Ml IKy^ T T 

iNUrll 1 

NU 
NX 
NPP 

NUMRUN 
OPT 


001 . 








2) State Variable Transition 
Matrix 












(read in by rows) 


F 




NX 


NX 


(10 F8.0) 


3) Control Variable Transition 
Matrix 

(read in by rows) 


G 




NX 


NU 


(10 F8.0) 


4) Exogenous Variable Transition 
Matrix 

(read in by rows) 


H 




NX 


NX 


(10 F8.0) 


5) Exogenous Variables for all 
Periods 

(read in by rows) 


Z 




NPP 


NX 


(10 F8.0) 


6) Initial State Variables 


X 






1 


(10 F8.0) 


7) Initial Control Variables 
for al 1 Periods 
(read in by rows) 


U 




NPP 


NU 


(10 F8.0) 


IL DATA FOR CALCV SUBROUTINE 


1) Vector of p.f. Weights. 


K 




9 


1 


(9 F8.0) 


2) Vector of p.f. Ratios 
(Note R(5) is wt for LPP 
X9,x-|o>xn) 


R 




8 


1 


(8 F8.0) 


3) Vector of Coefficients for 
ASF Needs 


C 




3 


1 


(3 F8.0) 


END CARD 




100 






(13) 
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DATA STRUCTURE FOR MULTIPLE RUNS 



I tern 


Var .Name 


Value 


# Cols. 


# Rows 


Format 


1 ) RUN card (as before) 

2) MAIN data (as before) 

3) CALCV data (as before) 












4) RUN card (as before except 
run number) 
e) 


NUMRUN 


>002 








5) Change code 

New F-Matrix follows 
New G-Matrix follows 
New H-Matrix follows 
New Z-Matrix follows 
New Initial -state Variables 

follow • 
New Ini tial -control Variables 

TO : 1 OW 

No (further) changes in MAIN 
data 

Following each change coded card, 

1*hp a nnpnnp'i "hp vyapiahlpc app 
uiic ajjprupr la uc vai lauic^ arc 

specified as before. 




1 7 

2 
3 
4 
5 
6 

7 

1 






^ i- 1 


6a) If there are no changes in 
CALCV data 

6b) If there are changes in CALCV 
data 

Complete set of CALCV data as 
before 


ICHNGE 
ICHNGE 


0 
2 






(n) 
(n) 


Repeat 4-6 
as needed 












7) END CARD 




100 






(13) 
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